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^ Abstract 

We present a Monte Carlo code dedicated to the simulation of bremsstrahlung arising in collisions of polarized electrons with thin 
CD target foils. The program consists of an electron transport algorithm taking into account elastic electron-nucleus scattering and 
,— ' inelastic collisions with target electrons as well as a treatment of polarized-electron bremsstrahlung emission. Good agreement is 
, found between the predictions of the electron transport code and data stemming from other simulation programs and experiments. 

f"| In addition, we present first results from the bremsstrahlung simulation which indicate a significant decrease in the degree of linear 

Q-i polarization of bremsstrahlung even for the thinnest gold targets considered. 

^ Keywords: Monte Carlo simulation, polarized electrons, bremsstrahlung polarization, target-thickness effects, hard x-ray 
O polarimetry 
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1. Introduction 

Bremsstrahlung arising from the interaction of an energetic 
electron with a screened nucleus potential, also referred to as 
ordinary bremsstrahlung, is one of the basic photon-matter pro- 
cesses and has attracted continuous interest both theoretically 
and experimentally during the last decades Of par- 

ticular importance is the study of the bremsstrahlung process 
resulting from polarized electrons as it reveals subtle details of 
the polarization transfer from charged particles to photons. The 
dependency of bremsstrahlung on the electron spin can be ex- 
pressed in terms of the so-called polarization correlations that 
have been systematically studied by Tseng and Pratt jkj] and 
were recently revisited in several theoretical works [QeI 01- 
Recent interest focussed on the polarization properties of the 
emitted bremsstrahlung photons with respect to the incoming 
electron polarization. In general, the use of polarized elec- 
trons should lead to a significant change in the degree of lin- 
ear polarization and a rotation of the polarization axis of the 
bremsstrahlung photons with respect to the unpolarized elec- 
tron case, where the photon polarization is solely determined by 
the kinematic and atomic parameters of the interaction. It is in- 
teresting to note, that a similar effect was recently predicted for 
the radiative recombination process in case of polarized highly- 
charged, heavy ions, see |Ht] for details. These theoretical stud- 
ies were mainly motivated by the development of novel Comp- 
ton polarimeters that enable efficient and precise measurements 
of linear polarization in the hard x-ray regime 
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During the last two years, such polarimeters were applied 
in a series of polarization-resolved bremsstrahlung measure- 
ments that have been carried out at the teststand of the polar- 
ized electron source SPIN iQj, 14] at the Technical Univer- 
sity of Darmstadt. In these experiments the polarization trans- 
fer was studied in collisions of polarized electrons with gold 
foils of about 100/ug/cm 2 thickness and at an impact energy 
of lOOkeV E El. This kind of investigations may open 
a route for 'complete' measurements where besides the en- 
ergy and angular information also the polarization properties 
of the involved particles are obtained [1]. However, a crucial 
point when drawing conclusions from the experimental data 
is the question to what extend the bremsstrahlung character- 
istics, in particular concerning the polarization, is altered by 
effects due to multiple collisions inside the target foil. While 
the influence of the target thickness on bremsstrahlung angular 
and spectral distribution was intensively investigated in several 
works H3, [3 El, S El Q, to the best of our knowledge no 
systematic study of target effects on the bremsstrahlung polar- 
ization was ever published. However, one can expect that the 
linear polarization shows a high sensitivity to target effects as 
the straggling of incident electrons followed by bremsstrahlung 
emission will lead to a superposition and, consequently, to a 
partial cancelation of different polarization orientations in the 
detector. 

In general, Monte Carlo simulations are well-suited for 
the determination of multiple-scattering effects on the proper- 
ties of processes where rigorous theoretical treatments assume 
single-collision conditions. However, in case of bremsstrahlung 
most of the general-purpose Monte Carlo programs available 
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(PENELOPE, Geant4, EGS5, etc.) use approximations that are 
limited to thick-target bremsstrahlung and/or ignore the effects 
of electron polarization on the bremsstrahlung properties . 

In this work we present the Monte Carlo code PEBSI 
(Polarized Electron Bremsstrahlung Simulator) that models the 
transport of polarized electrons through the target material in 
order to estimate the bremsstrahlung properties resulting from 
the finite thickness of the target foils. Here, the bremsstrahlung 
theoretical data are taken from fully relativistic calculations tak- 
ing into account electron and photon polarization. The devel- 
opment of the code was originally motivated by bremsstrahlung 
polarimetry measurements as described above. However, it 
might be applicable to a variety of studies dealing with polar- 
ized electron transport and bremsstrahlung emission. 



2. Polarized electron transport model 

In the energy region between a few keV up to a few MeV 
the electron transport in solid state targets is governed by two 
processes, namely the elastic electron-nucleus scattering and 
the inelastic scattering on bound target electrons, with the lat- 
ter leading to a successive energy loss of the incident elec- 
trons. In the case of high-Z targets, the Z 2 scaling of the scat- 
tering cross sections leads to a dominance of the elastic scat- 
tering process. Thus, for our application the elastic scatter- 
ing cross section should be treated as accurately as possible, 
while for the inelastic scattering even a rather approximate han- 
dling will not significantly decrease the validity of the electron 
transport model J23 , 24 , 25 ] . Note that we do not consider the 
bremsstrahlung process within the electron transport code as its 
contribution is negligible in the electron energy region of inter- 
est. Consequently, the mean free path A between two interac- 
tions in the target is given by 
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where p is the target density and cr e \ and cr inel denote the 
elastic and inelastic interaction cross sections per atom, respec- 
tively. The incident electrons are followed as long as they do not 
leave the target foil (transmission or backscattering) and their 
kinetic energy stays above a certain cut-off value E cut . 

The general structure of our Monte Carlo program is adopted 
from the approach presented in 1I26II . In the following we briefly 
discuss the implementation of both processes from Eq. [T]in the 
PEBSI code. In addition, results of the electron transport code 
are compared to data from experiments and a different Monte 
Carlo program. 

2.1. Elastic electron-nucleus scattering 

Here, we summarize the main aspects of the elastic electron- 
nucleus scattering process as it is implemented in the PEBSI 
code. For a much more detailed discussion we refer the reader 
to & 

In case of unpolarized electrons the differential cross section 
is given by the Mott equation 1 28] : 
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where 6 denotes the polar scattering angle and / is the am- 
plitude of the scattering wave with spin direction remaining un- 
changed while g is the amplitude of the spin-flip scattering wave 
as a result of the spin-orbit coupling. 

When polarized electrons are considered an additional de- 
pendence on the azimuthal scattering angle <p is introduced: 

£.(£) i(1+sm ,..,. 

where P denotes the orientation of the incident electron spin 
and h is the unit vector perpendicular to the scattering plane, 
which is defined by the directions of the incoming and the out- 
going electron. S(0) is called the Sherman function and is 
given by 
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with * indicating the complex conjugate. According to Eq.[3] 
incident electrons with a spin orientation having a transversal 
component with respect to the electron momentum will exhibit 
an anisotropic azimuthal scatter distribution with the degree of 
this asymmetry being determined by the value of S(8). This 
effect is used in Mott polarimetry where the degree of electron 
spin polarization is obtained by measuring the elastically scat- 
tered electron intensity distribution, see H29I1 for details. 

After the scattering process took place the new orientation of 
the electron spin P' is given by 
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As for the treatment of elastic scattering in PEBSI, the 
(s§)o va l ues ^ taken from the NIST database of electron 
elastic-scattering cross sections (version 3.1) 113 Oil while the S , 
T and U functions are derived from the Dirac equation using 
the Thomas-Fermi-Dirac atomic potential as described in [31]. 



2.2. Inelastic electron-electron scattering 

To model the inelastic scattering and also the energy loss of 
incident electrons we applied the commonly used approach of 
a hybrid model in which inelastic interactions are considered as 
to be either 'hard' or 'soft'. Hard interactions imply a signifi- 
cant amount of energy loss AE > Ethresh an d angular deflection 
sin 2 AO - AE/E, and are therefore explicitly treated as discrete 
collision events according to Eq. Q] On the other hand soft col- 
lisions are considered as a continuous process with the energy 
loss being proportional to the distance traveled by the electron 
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between two discrete interactions while the momentum and the 
spin orientation of the electrons remain unchanged. In this work 
we set Ethresh = 1 keV. The stopping power due to soft interac- 
tions is then obtained from 
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where is the CSDA (continuous slowing down 

V ds / CSDA v ° 

approximation) stopping power provided by the ESTAR pro- 
gram JH from NIST and (^f ) hard is calculated as 
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with the incident electron kinetic energy E and the fractional 
energy loss e = AE/E. The minimum and the maximum energy 
loss for a specific interaction are given by e m i„ and e max , re- 
spectively. The lower threshold e m i n is defined as the maximum 
value of two parameters, namely the energy that distinguishes 
hard from soft collisions ^thresh an d the binding energy of a spe- 
cific target electron, while the upper threshold e max equals either 
0.5 (free-free collisions) or 1 (free-bound). The effective cross 
section o"i ne i is the sum of the individual inelastic interaction 
cross sections being multiplied by the number of respective tar- 
get electrons as described in the following. 

If the binding energy of the target electrons can be ne- 
glected the electron-electron scattering process is described by 
the M0ller cross section for free electrons: 
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Here r e denotes the classical electrons radius, y is the rel- 
ativistic factor of the incident electron and nf ree is the number 
of target electrons considered as quasi-free. This approximation 
was applied to all electrons with binding energies /, < £Vesh, re- 
sulting in e m i n = E t hr es h/E. Note that besides the fact that Eq.[8] 
is not applicable for collisions where the energy transfer is in 
the order of the binding energy, the lower threshold e m ; n > in 
Eq. [7] is also necessary because the M0ller cross section leads 
to a divergence of the integral when e m i„ -»■ 0. We have to stress 
that because of this limitations the PEBSI code is not able to re- 
produce subtle features that are connected to low-energetic in- 
teractions and can only be treated adequately when using much 
more sophisticated models, see, e.g., B33I1 and II34I1 . However, as 
the energy resolution of the x-ray polarimeters mentioned above 
is in the order of 2 keV such features will not significantly alter 
the properties of the detected bremsstrahlung. 

For target electron shells i with binding energies /, > Ethresh 
we apply the Gryzinski cross section fl35ll : 
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Figure 1 : Inelastic cross sections as a function of electron impact energy mul- 
tiplied by the number of target electrons for collisions with gold atoms. The 
PEBSI cross sections are compared to various experimental and theoretical 
data: DWBA calculations - • ♦ * |3§], experiment - O EjJ O j4l|] □ B3l 
# El A m v Q. 
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with e, = Ij/E and the prefactor C2 given by 
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The R denotes a correction factor accounting for relativis- 
tic incident electron energies and/or relativistic binding ener- 
gies which was also introduced by Gryzinski. Here, we use a 
slightly modified R that was proposed in ll36ll . Though rely- 
ing on a purely classical collision model, the electron impact 
ionization cross sections yielded by Eq. [9] are often in reason- 
able agreement with data from experiments and more sophis- 
ticated calculations. However, in the near-threshold region the 
Gryzinski model tends to significantly underestimate the total 
cross section values, see 13711 . Therefore we replaced the total 
electron impact ionization cross section for the K-and the L- 
shell electrons by the semi-empirical formula given by Haque 
et al. Ik36l, l38ll . while keeping the energy differential behavior 
given by Eq.[9j 

Fig. Q] shows the resulting cross sections as a function of the 
incident electron energy. The PEBSI results are compared to 
data from DWBA calculations and several experiments. In gen- 
eral, reasonable agreement is found. Although the collisions 
with inner-shell electrons give a minor contribution to the elec- 
tron transport as it is dominated by the elastic scattering pro- 
cess, the cross sections can be used to model the resulting fluo- 
rescence x-rays which may alter the bremsstrahlung radiation. 
In the future, we also might take into account the secondary 
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electrons generated by electron-electron collisions in order to 
investigate their contribution to the bremsstrahlung. 

Note that in case of electron-electron collisions, we do not 
consider the electron spin polarization as the scattered electron 
distribution is altered only in the case where both collision part- 
ners are polarized. However, the incident electron spin orien- 
tation is altered during the collision leading to a successive de- 
polarization of the electron beam. At the moment this effect 
is taken into account by setting the electron to be unpolarized 
after the first hard collision. Naturally, this assumption overes- 
timates the effect of depolarization and we plan to implement 
are more realistic treatment in the near future. 

2.3. Verification of the model 

In order test the reliability of our electron transport model 
we compare the PEBSI results with predictions by a different 
Monte Carlo code as well as with experimental data. 



100 




10000 



Foil thickness [pg/cm ] 



Figure 2: Backscattering and transmission probability for electrons impinging 
on gold foils of various thicknesses. The PEBSI results are compared to pre- 
dictions by the PENELOPE program 1 45 ] . The lines connecting the data points 
are drawn to guide the eye. 

In Fig. [2] the percentage of backscattered and transmitted 
electrons is shown as a function of target thickness for dif- 
ferent incident electron energies. We simulated an unpolar- 
ized electron beam impinging at normal incidence on gold 
foils. The data from PEBSI are compared to predictions by 



the PENELOPE code (version 2008) ||45|], which is a general- 
purpose Monte Carlo program package for electron and pho- 
ton transport and is known to provide a very reliable trans- 
port algorithm for unpolarized electrons Il46ll . In order to per- 
form the comparison under almost identical conditions, both 
codes were adjusted to discard electrons with energies falling 
below £c Ut = 1 keV and the minimum energy loss treated as 
a hard inelastic collision was set to 1 keV (PEBSI) or to 1 % 
of the electron energy (PENELOPE), respectively. Moreover, 
PENELOPE was set to treat the elastic scattering process as 
detailed as possible (by setting the maximum angular deflec- 
tion that is treated with a multiple-scattering approach to zero). 
The generation of secondary electrons as well as the emission 
of bremsstrahlung was not taken into account. As seen, both 
programs yield nearly identical results. 
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Figure 3: Anisotropy parameter S e g of the scattered electron distribution in 
case of lOOkeV electrons being transversely polarized with respect to the elec- 
tron momentum. The PEBSI results for gold targets of various thicknesses are 
compared to experimental data 14711 . 



In case of electrons being transversely polarized with respect 
to the electron momentum, the azimuthal asymmetry of the 
elastically scattered electron distribution is characterized by the 
effective Sherman function S e ff(8). For infinite-thin targets one 
expects that S e s equals S and, consequently, the scatter distri- 
bution is given by Eq. [3]as discussed above. In Fig. [3j PEBSI 
results obtained for gold targets of various thicknesses are com- 
pared to experimental data which were measured for an electron 
energy of 100 keV [47]. The theoretical Sherman function for 
single-collision conditions is shown in addition. Although one 
finds a qualitative agreement between the data from simulation 
and experiment, a quantitative analysis would require a much 
more refined simulation taking into account the complete ex- 
perimental setup. For example, the decrease of the experimen- 
tal S e ff in case of the thinnest target and at backward angles 
might be due to electrons being backscattered from the beam 
dump behind the target foil and the dip near 90° was probably 
caused by the target holder. However, such a detailed study is 
beyond the scope of the present work. 



4 



3. Simulation of bremsstrahlung linear polarization 

As discussed above, the emission of bremsstrahlung radia- 
tion is not implemented as a part of the electron transport code. 
Instead, we treat the electrons like they were permanently emit- 
ting bremsstrahlung photons while we ignore the effect on the 
electron properties. The procedure is as follows: At every dis- 
crete interaction point according to Eq. Q] we calculate the ob- 
servation angle of a virtual x-ray detector with respect to the 
actual electron momentum. In addition, we estimate the elec- 
tron spin orientation with respect to the reaction plane. This 
information together with the electron energy is used to obtain 
the cross section and polarization properties of bremsstrahlung 
which was emitted along the straight path from the previous in- 
teraction point to the actual one and was heading in the direction 
of the x-ray detector. As the bremsstrahlung emission proba- 
bility is proportional to the distance traveled, the path length 
between both interactions is used as weighting factor. The 
theoretical electron-nucleus bremsstrahlung data is taken from 
tables which were calculated within a fully relativistic treat- 

r— i 

ment, see [5] for details. Note that we do not take into account 
the contribution of electron-electron bremsstrahlung, which is 
roughly a factor Z smaller than the ordinary bremsstrahlung. 
Scattering and absorption processes of the bremsstrahlung pho- 
tons inside the target are also neglected. The superposition of 
the radiation stemming from a large number of electrons along 
their complete tracks inside the target material then transforms 
the theoretical data for the single-collision case to effective 
bremsstrahlung properties taking into account the target effects 
as well as the position and solid-angle coverage of the virtual 
detectors. 
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Figure 4: The degree of linear polarization of the high energy end of 
bremsstrahlung plotted as a function of the observation angle and various target 
thicknesses. The PEBSI results for unpolarized electrons impinging on gold 
targets are compared to the underlying theoretical predictions for the single- 
collision case | 

As a first result, we present in Fig. |4]the degree linear polar- 
ization of bremsstrahlung stemming from unpolarized electrons 
impinging on gold foils for various target thicknesses in com- 



parison to the underlying theoretical predictions for the single- 
collision case. The PEBSI results are shown for the high-energy 
end (hu\E = 0.99), also referred to as the short-wavelength 
limit, of the bremsstrahlung distribution where the polarization 
is most pronounced. In order to visualize target effects only, 
the virtual detectors as well as the electron beam radius were 
assumed to be point-like. A clear decrease in the degree of 
linear polarization is found even for the thinnest target foil con- 
sidered. This finding is remarkable as for high-Z targets in the 
literature, thicknesses in the order of 50 yug/cm 2 are regarded as 
thin enough in order to minimize the effect of multiple colli- 
sions on the emitted bremsstrahlung radiation [3J. As a con- 
sequence, one can in general expect significant target effects 
on the bremsstrahlung linear polarization for self-supporting 
high-Z target foils that are typically used in experiments. 
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Figure 5: Bremsstrahlung linear polarization as a function of the photon energy 
for an observation angle of 120° with respect to the direction of the incident 
electron beam, see text for details. 



In fig. [5] the degree of linear polarization is shown as a func- 
tion of the bremsstrahlung photon energy for an observation an- 
gle of 120°. For target thicknesses above 100 fig/cm 2 , the depo- 
larization due to target effects decreases with increasing photon 
energy at the high-energy end of the bremsstrahlung distribu- 
tion. This feature is caused by the fact that electrons with large 
deflection angles have such long tracks inside the target that 
their energy loss becomes significant and, consequently, after 
a certain track length these electrons have not enough kinetic 
energy to contribute to the highest bremsstrahlung energies. 
However, they can still emit bremsstrahlung with slightly lower 
photon energies and here their large deflection with respect to 
the incident electron direction leads to a strong depolarization 
when summing over the radiation emitted by a large number of 
individual electrons. This in contrast to thinner targets where 
only the small fraction of electrons with deflection angles very 
close to 90° (when assuming normal incidence of the incom- 
ing electrons) is likely to stay long enough inside the target to 
suffer a significant energy loss. From this results one can ex- 
pect an almost complete depolarization of bremsstrahlung when 
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high-Z targets with thicknesses in the order of 1 fjm (equals 
1932/ig/cm 2 in the case of gold) are used and photon energies 
at least a few keV lower than the short-wavelength limit are 
observed. 

As for the case of bremsstrahlung arising from polarized 
electrons, the development of the code is still ongoing. How- 
ever, preliminary results already indicate that the PEBSI code 
is able to qualitatively reproduce the bremsstrahlung properties 
that were obtained in recent experiments. 
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